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Abstract 

Asymptotic formulae for the mechanical and electric fields in a piezoelectric body with a small 
void are derived and justified. Such results are new and useful for applications in the field of design 
of smart materials. In this way the topological derivatives of shape functionals are obtained for 
piezoelectricity. The asymptotic formulae are given in terms of the so-called polarization tensors 
(matrices) which are determined by the integral characteristics of voids. The distinguished feature 
of the piezoelectricity boundary value problems under considerations is the absence of positive defi- 
niteness of an difi'erential operator which is non self-adjoint. Two specific Gibbs' functionals of the 
problem are defined by the energy and the electric enthalpy. The topological derivatives are defined 
in different manners for each of the governing functionals. Actually, the topological derivative of the 
enthalpy functional is local i.e., defined by the pointwise values of the governing fields, in contrary to 
the energy functional and some other suitable shape functionals which admit non-local topological 
derivatives, i.e., depending on the whole problem data. An example with the weak interaction be- 
tween mechanical and electric fields provides the explicit asymptotic expansions and can be directly 
used in numerical procedures of optimal design for smart materials. 

Key words: Piezoelectricity, polarization matrix, asymptotic analysis, electric enthalpy, topolog- 
ical derivative, optimum design, shape optimization. 
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1 Introduction 



The paper is devoted to the asymptotic analysis of boundary value problems for coupled models. The 
coupling occurs between the mechanical part which takes the form of the linearized elasticity and governs 



the stress-strain state of the body, and the electrical part which describes the electromagnetic field in 
the body. 

From the view point of applications, piezoelectric materials are of common use in electromechanical 
sensors and actuators, e.g., ultrasound transducers in medical imaging and therapy, force and acceleration 
sensors, positioning sensors, surface acoustic wave filters, still with the growing range of applications in 
modern technology. Their mode of action is based on the piezoelectric effect, that couples the electrical 
and mechanical behavior of such materials. For the optimal design of piezoelectric devices, efficient 
numerical procedures for shape and topology optimization should be still developed. In the modern 
theory of shape optimization it is required that the derivation of shape and topological derivatives of 
shape functionals to be optimized is performed beforehand. From one side, the derivation of shape 
gradients of integral functionals in smooth domains [49] and non-smooth domains [28] (cf. [lQl[4T]) has 
become a standard procedure. There is no major difficulty to perform such a shape sensitivity analysis 
for the elliptic boundary value problem under considerations. However, the boundary value problem 
in piezoelectricity cannot be posed in such a way that it simultaneously is formally self-adjoint and 
possesses a semi-bounded quadratic form. This specific feature makes the problem more involved from 
the asymptotic analysis point of view compared to the pure elasticity or pure electricity boundary value 
problems. In addition, the general case of inhomogeneous and anisotropic body is considered, which also 
requires for additional and new technicalities in asymptotic procedures which is the main subject of the 
paper. In particular, different formulations of the piezoelectricity problem (cf. Sections 12. 3[ 12. 41 13. 3p lead 
to two definitions of the polarization matrices which differ one from another by its properties. Moreover, 
only the electrical enthalpy, which is but the governing functional for the piezoelectric media (see, e.g., 
[12] [T5l |46] ) admits the topological derivative dependent on local characteristics of mechanical and 
electrical fields. Other shape functionals, especially the energy functional, get the topological derivatives 
dependent on the global characteristics of mechanical and electrical fields. This acquired trait raises the 
natural question on the properties of material derivatives for piezoelectricity in the framework of the 
shape sensitivity analysis with smooth or non-smooth boundary variations, it is clear that the result 
could be of the same nature, since the topological derivatives can be identified from the first order shape 
gradients by a limit passage e.g. in elasticity, |5j (cf. also [T4]). 

In the paper, we restrict ourselves to the asymptotic procedures of singular domain perturbations 
which allow us to obtain, in a natural way, the topological derivatives of shape functionals. In principle, 
the method developed here can be generalized to characterize the influence on solutions of the non- 
smooth boundary variations, therefore, we can derive the shape gradients even in such a case, e.g., for 
small defects located close- by the boundary, including micro-cracks (see [42]). 

Without entering into details, but with the strong practical implications in mind, we can claim 
that some possible applications of shape optimization in the field concern the design of electro-acoustic 
transducers which are constructed with piezoelectric actuator-patches and capacitative micro-machined 
ultrasound transducers. The task for optimal design for a class of electrostatic-mechanical-acoustic 
transducers can be e.g., the topology of electro-acoustic material and the topology of the electrode- 
layers, in order to achieve a maximal acoustic pressure, or a maximal acoustic energy in a specific 
sub-domains of the hold-all-domain. We refer the reader e.g. to [ISKIIKIS] for modeling of piezoelectric 
materials, to [18] for material tensor identification for such materials, and to for control issues. 

Our aim is a possible application in shape optimization, thus we introduce the so-called topological 
derivatives of shape functionals for piezoelectric materials. It seems that the models are not up to now 
used in applied mathematics for the purposes of shape optimization, although the smart materials are 
of common use in the engineering practice. In shape optimization, the modern approach to numerical 
solution, requires the preliminary knowledge of explicit formulae for shape gradients [49] as well as of the 
topological derivatives [50l HO] [TOj . These formulae are used in the level-set-type methods which model 
the geometrical domain evolution by a zero-level set of solutions to non-linear hyperbolic equations of 
the Hamilton- Jacoby type. The shape gradient are present as the coefficients of the equations, and the 
topological derivatives are used to improve the values of the shape functional under consideration by 
the appropriate topology changes, e.g., for the minimization of the shape functional, the minima of the 
topological derivative of the functional indicate the location of a new hole in the geometrical domain 
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2 Preliminaries. Problem formulation and description of re- 
sults 

2.1 Shape optimization in piezoelectricity 

This paper is motivated by the fact that, among numerous pubhcations on shape optimization, shape sen- 
sitivity analysis for piezoelectric bodies does not exist, although piezoelectric materials are of extremely 
wide usage in the modern technologies, one can think of a simple lighter, available in any supermarket, or 
an elaborated computer work-station in a university. One, and definitely not the only one, distinguishing 
feature of such smart materials implies an easy energy transfer in both directions from mechanical fields 
to electric fields. The mathematical modeling of such a phenomenon leads to serious complications of 
analysis for governing partial differential equations because the corresponding boundary value problem 
is not formally self-adjoint in contrast to the boundary value problems for purely elastic bodies or purely 
electromagnetic media. This fact requires for the development of new mathematical tools and a careful 
choice of the cardinal shape functional while neglecting of non-self-adjointness provokes mistakes in both, 
mathematical formulae and physical interpretation of the obtained results (see Remark 1211 below) . 

Introduced in [500, the topological derivative T(m°; wi) of a shape functional J is intended to describe 
the change of the functional J due to nucleation of holes or voids and allows to extend possible variations 
of the shape in an optimization process [2l[8l[9] in comparison with classical tools (cf. [49l [5l [6]). 

J(m''; ^{h)) = J{u; Vl) + h'^rixL- cc^i) + o(/i''), h ^ +0, (2.1) 

In (|2.ip . ft, > is a small parameter, i.e., the diameter of the opening cj^ in the entire body C M", 
v}^ and u are solutions of the boundary value problem in r2(/i) = \ cJ/i and fi, respectively, and the 
exponent k > depends on the space dimension n and boundary conditions imposed on the interior 
(n — l)-dimensional surface dujh- 

Asymptotic analysis of elliptic problems in singularly perturbed domains, e.g., methods of matched 
and compound asymptotic expansions (cf . [16j and |28| , respectively) , has become the most appropriate 
and relevant to obtain almost explicit formulae for the topological derivatives as it has been demonstrated 
in [40l|4T] and others. We also mention books [32l|3] where the subject is studied, to some extend, from 
physical and numerical point of view. 

Strangely enough, only self-adjoint problems were heretofore examined carefully, although the full- 
blown approach in [28| can work for general boundary value problems for elliptic systems. In this paper 
we partly fill this gap by adapting formula (|2.ip to certain shape functionals for a piezoelectric body. 

The piezoelectricity problems admits two different formulations with non-symmetric and symmetric 
but non-semibounded quadratic forms, the energy and electric enthalpy functionals, respectively. By 
means of the Lax-Milgram lemma, the first formulation furnishes the existence and uniqueness result. 
At the same time, the topological derivative of the energy functional is a non-local characteristics of the 
piezoelectricity solutions in contrast to the pure elasticity problem (see Remark below) , while the 
absence of this intrinsic property is not caused by an incorrect definition (|2.ip but occurs as well for 
the energy release rate in mechanics of cracks for piezoelectric media (sec Remark [21] again) . The fair 
explanation, we refer the reader to [S2] for the modeling issues, of the latter refers to the electric enthalpy 
as one of Gibbs' functional obtained from the energy functional by the partial Lagrange transform on the 
electric components. This is the electric enthalpy £(it''; ri(/i)) (see the definition in (|2.20p ). that governs 
the mechanical electric state of the piezoelectric body J7(/i) and, therefore, the second formulation 
becomes variational and provides the clear interpretation of the topological derivative 7£(u; wi) in 

£{uh;n{h)) = £{u;n) + h^T£{u;LUi) +0{h^), (2.2) 

^Actually, asymptotic formulae of type l|2.1|l together with the whole asymptotic scries for energy functionals under 
various singular boundary perturbations had been derived much earlier in |27| . although the notion of the topological 
derivative is due to I50| . 
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The indicated peculiarity of the piezoelectricity problem crucially influences topological derivatives 
of other shape functionals, too. For example, the traditional adjoint state (cf. OEHllSn]) has to be 
found out in the formally adjoint boundary value problem that occasionally underlines its name. 

All the above observations lifts the piezoelectricity problem on the top of the list of unsolved problems 
in shape optimization, it seems that even the classical formulae for material derivatives, which are not 
under consideration in the paper, ought to be revisited. 

2.2 Methods of asymptotic analysis 

Nowadays there exist several methods to construct asymptotic expansions of solutions to elliptic bound- 
ary value problems in domains with singular perturbations of boundaries. First of all, we mention two 
methods, namely, the method of mathed asymptotic expansions and the method of compound asymp- 
totic expansions (cf. monographs [16] and 01] j respectively), which in general appear to be of the same 
power. Indeed, based on different asymptotic procedures, the matching procedure and the procedure of 
of rearrangement of discrepancies, they result in asymptotic expansions which differ at the first sight 
one from another, but can be readily transformed one into another (we refer to the introductory chapter 
2 in [28]). By the way, we silently use this transformation while presenting at the end of Section [3751 
an alternative way of presentation the asymptotic form of the derived solution. The method of com- 
pound asymptotics is employed throughout the paper for two reasons. First, the results given in [28] 
are obtained in relatively general formulation which includes the elliptic systems of partial differential 
equations not necessarily formally self-adjoint (cf. discussion in the preceding Section [^TT|) . On the other 
hand, in the monograph [16] the results are established exclusively for the scalar second-order elliptic 
equations in the divergence form. Second, the method of compound asymptotic expansions is carefully 
matched with theory of elliptic problems in domains with conical outlets to infinity, specifically the ex- 
terior domains (cf. [(3^; Ch.6], and [35]) while for our purposes this theory is used further to introduce 
and investigate the polarization matrices in piezoelectricity. 

Since the problem under studies is geometrically very specific, i.e., it concerns only one small opening 
inside of a domain in M^, the other methods of asymptotic analysis can be employed. In Remark [TSl 
(2) we mention the case of piecewise constant coefficients which makes suitable an asymptotic analysis 
of the equivalent boundary integral equations obtained from fundamental solution (cf. [3] with similar 
results in elasticity), although no fundamental matrix is known in piezoelectricity. The other possibilities 
include among others the homogeneization technique relying on the so-called delute limit (see [T7] , [3T] , 
|26| and many others). However, in our opinion, the method of compound asymptotic expansions is 
still the most appropriate tool in piezoelectricity in order to investigate asymptotic properties of shape 
functionals. 

2.3 Constitutive relations in piezoelectricity 

Let C be a piezoelectric body with the Lipschitz boundary dft and the compact closure = flUdfl. 
Using the matrix/column notation (cf. [Ml[33]), we regard the displacement vector as the column 
~ {u^,u^,u^)^ where is the projection of u on the Xj-axis of the fixed Cartesian coordinates 
system x = {xi,X2,X3)'^ and T stands for transposition. Together with the electric potential u^, the 
displacements compose the column u = {u^ , , , u^)^ of height 4. The strain column 

e^«^) = {sr,,s^„e^,,V2e^„V2e'^„V2er,V (2-3) 

consists of the Cartesian components e'^^, = ^{9jU^ + dkU^) of the strain tensor and takes the form 
e'^{u'^) = D'^{V^)u'^ where 

( di 
V 





2-1/2^3 

^3 2-1/2^2 





(2.4) 
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We introduce the column e{u) = {e^{u^)^ ,e^{u^)^)^ where e^(u^) = V^it^ is the electric strain 
column, taken with the sign minus, and D(Vx) implies a (9 x 4)-matrix of the first-order differential 
operators, 



eiu) = D{V,)u, D{V,)' = [ Q^^^^ ), 0== (0,0,0). (2.5) 



We also assemble the column a{u) of height 9 from the stress column (t^{u^) of structure (|2.3p 
and the electric induction column cr^(w^) = (cr^ , erf , erf )^. In this manner, the constitutive relations of 
piezoelectricity (see [l2 l [T5l l46] ) 

^ ^MM^M _ ^ME^E^ ^E ^ ^EM^M ^ ^EE^E (3.6) 

can be rewritten as follows: 



cr(u) = Ae{u), (2.7) 



where the matrix A of size 9x9, 



(2.8) 



is formed by the symmetric and positive definite (6 x 6)- and (3 x 3)-matrices and A^^, respectively 
the elastic stiffness matrix and the dielectric permeability matrix, and the blocks ^4^^ = (A^^)'^ of 
piezoelectric moduli. We emphasize that, by its physical nature, the matrix (|2.8p is not symmetric 
provided the (6 x 3)-block A^^ is not null, i.e., the mechanical and electric fields interact. 
The state of the piezoelectric body il is described by the mixed boundary value problem 

D{~V,yAix)DiV,)uix) = fix), X e (2.9) 
D{n{x)y A{x)D{V,j:)u{x) = g{x), x € T^, (2.10) 

u{x)^o, xer.u = dn\T,, (2.11) 

where n = {ni , ?i2 , '^a)^ is the unit vector (column) of the outward normal. On the right hand-side of the 
equations (|2.9[) and (|2.1G[) . we have the volume forces /'^ = (/i", /2''j /s")^ and the surface mechanical 
loading = {gi^ , , g^^)^ together with the volume and surface g^ electric charges. The Dirichlet 
conditions ()2.11|) mean that the body is mechanically clamped over the surface and in contact with 
an electric conductor. As usually, = and, if the surface is in contact with a dielectric medium, 
i.e., vacuum, we are to put gf = 0. 



2.4 Solvability of boundary value problem 

Let us assume that mes2r„ > and / G L^(ri)'*, g E i^(ro-)^ where i^(S) denote the Lebesque space 
with the intrinsic inner product ( , )h and the superscript 4 indicates the number of components in the 
vector functions / and g. Notice that the subscript is always omitted in our notation for inner products 
and norms. 

The integral identity (cf. [53]) serving for problem (|2.9|) - (|2.1ip . reads as follows: 

Qiu,v;n) := {AD{V ,)u, D{V ,)v)n = if,v)n + ig,v)r^, v e H\n;ru)\ (2.12) 

where H^{^; r„) denotes the Sobolev space of hmctions vanishing at Tu- The left-hand side of (|2.12p is 
understood properly provided entries of the matrix A are measurable and uniformly bounded functions 
in Q. In addition, for almost all x € n, we assume the symmetry and positivity properties 

A^^ix) = A'^^ix)^, A^^ix) = A^^ixy, A^^{x) ^ A^^{x)'^ , 

CM\a^\^ < {a'^yA^^ix) < CmIo^P, € M^^ (2.13) 

ceIo^P < {a^^A^^ix) < CeIq^P, e R\ 
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where cm, Cm and C£, Ce are positive constants. We emphasize that no positivity restriction is imposed 
on the piezoelectric moduU in A^^. 

Ahhough in the case A^^ ^ the sesquihnear form Q(-,-;ri) cannot be an inner product on the 
Hilbert space tl^i^^Yuf^ due to the wrong sign on A^^ in (|2.8p . the Lax-Milgram lemma ensures the 
following assertion because of the formula 

Qiu, u- n) = (A^^Z?^(V,)u^, i?^(V,)7.^)jj + (A^^V.u^ V,uE)n > c||u; H\n)\\^ (2.14) 

caused by the Poincarc inequality for and the Korn inequality for (see [ZlEl] and others). 

Proposition 1 Under the conditions (|2.13p . (j2.14p . the problem (|2.12p admits a unique solution u £ 
H^{^l^,Tu)'^ , and the following estimate is valid: 

\\u;H\n)\\<cn{\\f;L\n)\\ + \\g;L\r,)\\). (2.15) 

Unfortunately, the problem ()2.12p is non variational. Indeed, the energy functional U, 

U{u; VL) = '^{AD{V^)u, D{V,)u)n - A{u- U), (2.16) 
A{u-n):^U,u)n + {9.u)T,, (2.17) 
is but the sum of the mechanical and electric energy functionals 

u'^iu^'-n) = i(A^^i?^(v,)^.^,i5^(v,)«'')n - (.r,«^)n - {g''.u'')r^, (2.18) 

U^{u^;n) = ]^{A^^V,u^,V.,u%, - {f,u^)n - {9^,u^)r^, (2.19) 

while a stationary point of (|2.16p becomes a solution of the problem (|2.12p with the block-diagonal (9 x 9)- 
matrix diag{A^^ , A^^} , i.e., the variational problem docs not accept an interaction of the mechanical 
and electric fields (cf. an example in Section IT^ . 

It is known (see, e.g., [HI]) that the electric enthalpy £, 

£{u; n) = i(^(_)i?(V,)?/, D{V,)u)n - n{u; Q), (2.20) 
n(u;n) = (r,^.^)o + (3^,ii^)r. - {f!,ul)n - {glu^r^, (2.21) 
gives rise to the variational formulation of the piezoelectricity problem 

Q^_)iu,v;n):^iA^_)D{V^)u,D(y,)v)n=niv;n), v e H\Q;T,)\ (2.22) 
where the matrix is composed from blocks in (j2.8p as follows 

^(-) - 1 ^EM _^EE I [Z.26) 

The matrix (|2.23p . in contrast to the matrix A, is symmetric, however, neither matrix (|2.23p . nor the 
quadratic form on the left-hand side of (|2.22p is positive definite. Thus, a solution u G H^{fl;Tcr)'^ is a 
stationary point of the functional (|2.20p but u cannot be any minimizer of the electric enthalpy £{u] f2). 

The integral identity (|2.22p with the test function = {v^,v^,v^,—v^) transforms into the 

problem (|2.12p . The inverse transformation is also available. These facts prove that the problem (|2.22p 
inherits the unique solvability from (|2.12p in Proposition [1] 

Remark 2 The integral identity is formally obtained by the multiplying system (|2.9p with v scalarly and 
integrating by parts. Using as the multiplier, one arrives at (j2.22p . This explains the equivalency 
of the problems. 
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Remark 3 In Section \3.3\ we shall outline a different way to modify the piezoelectricity problem in order 
to study properties of the mechanical and electric fields on the base of known results. This approach 
is related to studies \31}/ . J^, and others, on variational formulations of elliptic problems describing 
processes in media with complex material coefficients (phase changes due to material properties). In 
fact we need this technique only for one reason, to maintain the so-called polynomial property \34\ \36^ . 
thus the results i31}/ . are not applied in the paper. We recall that the polynomial property allows 
to describe all required attributes of the exterior boundary value problem of piezoelectricity by simple 
algebraic calculations (cf. review \35^). 

The electric enthalpy is but the difference of elastic energy (|2.18|) and electric energy (|2.19p . Expres- 
sion (|2.17p implies the external work. Being the difference of the mechanical and electric external works, 
the component TZ{u; Jl) of the electric enthalpy has no physical meaning as a whole. Nevertheless, in 
Section 14721 we shall observe that asymptotic formulae for £{u] O) become meaningful while the analogous 
formulae for U(u;Q) look rather queer. 

2.5 Structure of the paper 

In Section [3] the asymptotic analysis of the piezoelectricity problem for the body n{h) with a small void 
TJh is performed (see p.ip ). The applied here asymptotic procedure [551 Ch.4] requires for introduction 
of an intrinsic integral characteristics of the void uJi in the homogeneous piezoelectric space M^, the 
polarization matrix M{A'^,uj) of size 9x9 (see formulae p.38p - p.40p ). Theorem [T2] establishes general 
properties of the polarization matrix, see also (|4.69p for the case of weak interaction between mechanical 
and electric fields. The polarization matrix appears in the asymptotic expansion of the boundary layer 
term at infinity that also permits in Section [3. 51 to complete the asymptotic ansatz of the solution to the 
piezoelectricity problem in fl(h). The asymptotics constructed in Section[3]is justified in Section [4. II In 
Section [4.2l the asymptotics of the energy and electric enthalpy functionals are analysed, while in Section 
14.31 rather arbitrary shape functional is considered and the corresponding adjoint state is detected. The 
paper is completed by inquiring into a piezoelectric body with a weak interaction of the mechanical and 
electric fields. All asymptotic formulae derived in the paper are made more explicit in such a case due 
to the fact that for pure electricity and pure elasticity the polarization matrices are known explicitly for 
many canonical shapes (see, respectively, [35], [SSKHIIS] and others). 

3 Asymptotic analysis 

3.1 The problem with an interior singular perturbation in the domain 

Let ijj be an open set in R""* with a Lipschitz boundary and a compact closure. We assume that both fl 
and uj contain the coordinate origin O. Given a small dimensionless parameter h G (0, Hq], we introduce 
the sets 

ujh = {a; : C := h^'^x G uj}, n{h) = 17 \ ZU,,. (3.1) 

The bound ft.o > is chosen such that oJ/i C 17 for h <E (0, /iq]. By rescaling, we reduce a characteristic 
size of n and lo to the unit and make the coordinates x and ^ dimensionless. 

Supposing Q{h) to be a connected set, we consider the piezoelectricity problem in the domain ^{h), 
namely, 

D{-V^)'^A{x)D{\/,)u\x) = fix), X e n{h), (3.2) 
Dinix)yA{x)D{y,)u'\x) = g(x), x e T,, (3.3) 
D{n''{x)y A{x)D{V^)u''{x) =0, xe duh, (3.4) 
u^\x) 0, a; e r„. (3.5) 
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In (j3.4p . n'* stands for the outward normal on duJh- Since the Neumann conditions arc imposed on the 
boundary of loh, there is no traction on dojh and the opening uTjl is fiUed with a dielectric medium. This 
problem, of course, ought to be reformulated as either integral identity (|2.12p . or (|2.22[) in the function 
space {VL{h)]T u)"^ , hence 

Q{u'\v'^-n{h)) = U,v%,^^) + {9.v^)r., G H'{n{h)-T^)\ (3.6) 

Proposition [1] remains valid for the problem (j3.6p in the domain Q,{h). 

For ft, = 0, the opening LJh disappears and the singularly perturbed problem (j3.2p - (j3.5p becomes 
the original problem (|2.9p - (|2.1ip . In order to describe the behavior of the solution E H^{^l{h);Tu)^ 
as ft ^ +0, we have to assume an additional smoothness of the matrix A, for example, in the ball 
B_R = {.T : |a;| < R} the inclusion 

is valid, where C'^'"(S) is the Holder space with the standard norm 

fc 

||z;;C^-"(S)|H5]sup|Vi«(a;)|+ sup |x - yr"|V^«(x) - Vjz;(y)| 

j=l ^^'^ x^y^i'B. 

and V^w denotes the family of all order k derivatives of v. Since the matrix differential operator 

L{x, V,) - D{-V :,)^ A{x)D{V ^) (3.8) 

is elliptic (see Section [3731 below) . a solution u E H^{Mi^)'^ of system (|2.9p in Mji with the right-hand side 

/eC°'"(B«)4, aG(l/2,l), (3 9) 

falls into the space C^'"(Bij')^ for any R' £ (0, i?). This fact is due to local estimates of solutions to 
elliptic systems [1]. Note that p.9p provides the estimate 

\fix)- fm<c\xr, x^^R- (3.10) 

We also need the Taylor formula 

\u{x) - d{x)a - D{xfe° - U{x)\ < c\x\^+", xGEr,, (3.11) 

where _D(a;)^ is the matrix in (|2.5p under the substitution V^; i-^ a;, 

e" = D{V^)u{0) e M^ (3.12) 

d{x)a with a G implies a rigid motion in the mechanical component and a constant potential in the 
electric one, 

1 -2-1/2x3 2-1/2^2 \ 

1 2-1/23.3 -2-i/2x-i . (3.13) 

1 -2-i/2a;2 2-^l^xi / 

We emphasize a similarity of the matrices D^{x)^ and d^{x). Finally, U in (|3.1ip is a quadratic term, 
i.e., 

U(tx) = t^U{x), t>0, xeR^. (3.14) 



8 



Remark 4 The factor \/2 is present in the strain column (|2.3|1 in order to equalize the natural norms for 
tensors of rank 2 with the norms of corresponding columns of height 6. As a result, an orthogonal trans- 
formation of the Cartesian coordinate system x implies the orthogonal transformations for all columns 
introduced to replace tensors (see, e.g., {23 Ch.2]). By the factor 2-^/2 in ({3^3)1 . we also achieve the 
velatioTis 

D{V^)D{xy = I9x9, D{V^)d{x) = ©9x7, ,o 1 

d(V,)Mx)U=0 = l7x7, rf(V,)^i?(x)^U=0-O7x9, ^ ' 

where Inxn o-nd Omxn stand for the unit and null matrices of size nxn and mxn, respectively. Notice 
that p.l2[) follows from the first couple of the relations (jS.lSp and our way to write the Taylor formula. 

By p.7p . we particularly obtain 

3 _ _ 

A{x) = A° + ^ XjA^ + A{x), \Apg{x)\ < c|.t|2, x e Mr, (3.16) 

with the constant (9 x 9)-matrices A^ so that matrix (|3.8p of differential operator gets the decomposition 

L(x,V,) =L"(V,)+i'(a-,V,)+Z(x-,V,). (3.17) 
Inserting the Taylor formula for u into the equation (|2.9p and using p.l6p yield 

3 

L"(V,)t/(x) - ^ D{e,)'^A^£" = /(O). (3.18) 
i=i 

Here cj = (Sj^i, 6j 2, Sj,^)^ . Since U is quadratic in a; (see p.l4p ). the first term on the left hand-side is 
independent of x. 

Remark 5 To guarantee formulae p.lOp and p. lip with a G (0,1/2), we could assume f G H'^{Mjf )^ 
while deriving u S iJ^(Bfl/)'^ from local estimates for solutions of elliptic systems (see ITI). This is due 
to the Sobolev embedding theorem H^^^ C C''" in R'^ for any a e (0, 1/2). However, in Theorem ] 19\ and 
Remark \20\ we shall see that we really need a > 1/2. The latter requires, for example, f G H^(Mii)'^, 
and, therefore, we prefer here to use the Holder scale. 



3.2 The asymptotic ansatz 

Based on general results in [28j on the asymptotic structure of solutions to elliptic boundary value 
problems in a domain with singular perturbations of the boundary, we accept the following asymptotic 
ansatz for the solution u'^ of problem p.2p - p.5p : 

u''{x) = u{x) + x{x){hw\^) + h^w^i^)) + h^u{x) + ... (3.19) 

Here u is a solution of the limit problem (|2.9p - (|2.1ip . and w'^ are terms of the boundary layer type, 
and u is the main regular corrector. The boundary layer terms are treated in Sections 13.31 and 13. 4|, 
respectively, and the regular corrector in Section [3751 below. The cut-off function x S Cj?°(J7) is equal 
to one in the ball B/;/3 and null outside M2R/3 so that, now, we fix Hq > such that ujh C for 
h€{0,ho]. 

Remark 6 The boundary layer solutions and uP' are constructed in Sections \3. 3\ and \3.4\ respectively, 
along with their decompositions at infinity. In Section \S.5\ the main terms of the decompositions, compose 
the right-hand side of a problem of type p.9p - (|2.1ip for the regular solution u. 
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In view of (j3.ip . the coordinate dilation x = h^^x removes the boundary dfl close to infinity 
and the formal limit passage h — > +0 makes the exterior domain S = M'^ \ ti; from the nucleated domain 
Moreover, the decomposition (|3.17[) yields 

L{x, V^) = L{h^, h~^V^) = h-^L"{V^) + h'^L'i^, Vj) + . . . (3.20) 

Similarly, for the Neumann boundary operator N^{x,Wx) on the left hand-side of (|3.4p . we have 

N'^ix, V,) = /i-iiV°(C, V^) + /i°7V'(e, V^) + . . . (3.21) 

where 

3 

ArO(e, V,) = D(n"(e))^A°D(V.), 7V'(C, V^) = D(n-(C))^ ^ e.^^'l>(Ve), (3.22) 

and is the unit vector of the outward normal on duj. 

Let us derive the exterior boundary value problems for and w'^. First, we insert the ansatz (j3.19p 
into p.2p . make use of the expansion p.2ip . and collect coefficients written in the fast variables ^ for 
similar powers of the small parameter h. As a result, we obtain systems of differential equations in S 
for and w'^ (see p.23p and (|3.24p below). Second, we calculate the discrepancy left by the leading 
asymptotic term u{x) in the boundary conditions (|3.4p . Namely, by means of (|3.1ip . (|3.2ip . we derive 
that 

N\x, Vx)u{x) = Din^iOV ( ^" + ^ E j + + ■■■ 



Finally, we write the problems 



LO(V5)t(;i(e)-0, eeS, 



(3.23) 



and 



(3.24) 



3.3 The exterior problem in piezoelectricity 

The polynomial property |34[|35| of a formally self-adjoint system of differential equations delivers plenty 
of results for the exterior boundary value problem in S such as the ellipticity, the solvability, asymptotic 
expansions of solutions, and intrinsic integral characteristics, i.e. the polarization matrices (see |391 
Ch.6], [3S1 [37] and [ID] in shape optimization). As it has been mentioned, the piezoelectricity system 
(|2.9p is not formally self-adjoint, however, introducing the imaginary potential mf (see [3SJ Example 
1.13]) and the column = {u^ , , , iu^)^ brings the sesquilinear form 

g(,)(u(,),«(,);S) = (A°)Z?(Ve)u(,),i?(Ve)i;(,))H (3.25) 

where i is the imaginary unit and A'^-^ stands for modified matrix (|2.23p . 



^ „ /10Em' ^OEE I — ^?fle) + *^am) I (3.26) 



^OMM j^OME 
HO - i^OEM'^ ^OEE J - ^(Re) -1- *^(/m)> 

while both ^(^g) and ^("/m) are real symmetric and ^(/jg) is positive definite. The sesquilinear form 
p.25p is not Hcrmitian in the case A^^^ ^ Oexs, but it enjoys the polynomial property [5H [5^ [55] : 

g(i)(u(,),U(i);T) = G "Plx, (3.27) 
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where T is any domain in M.^ and V = {p : p{x) — d{x)a, a G C^} is a polynomial subspacc of dimension 
7 generated by the matrix in (|3.13p . 

The above observations made in [36l [35] and the investigation scheme |39l Cli.6] provide all results 
we formulate below with exception for the polarization matrix and here the most attention is paid to 
this integral characteristics of the opening uJ in the homogeneous piezoelectric space. 

Let VqCE.) be the Kondratiev space [5D] obtained by the completion of the linear space C^(S) 
(infinitely diffcrcntiable functions with compact supports) with respect to the Dirichlet integral norm 
II V^w; L^(S)||. Applying the one-dimcnsional Hardy inequality in the radial variable p = |^|, we use the 
equivalent norm 

\\w; V,\E)\\ i\\V,w; L^m'' + Wp''^; L^m')'^'- (3.28) 

The problem (|3.23|) with the right-hand side g G L'^{duj)'^ in the Neumann boundary conditions can 
be reformulated as the integral identity, similarly to (|2.12p 

iA^D{V^)w, D{V^)v)s - (g, v)a^, v e Vo\E)\ (3.29) 

Proposition 7 For any g £ L'^{doj)'^, the problem (|3.29[) has a unique solution w G VqCE)^ and the 
estimate ||ti;; Vq'"(S)|| < cHg; L^(9w)|| is valid. 

Although duj and g are not smooth, the solution w in Proposition [7] is infinitely differentiablc outside 
of any neighborhood V of the set tU (recall the local estimates in [1] mentioned above) . To describe the 
behavior of w(^) as p oo, we introduce the fundamental matrix ^(x) of size 4x4 for the operator 
L^CVj) in R'^ (see [UllIS]). This matrix is positive homogeneous of degree —1, namely, 

$(tO = t-^'^i^), t > 0, ^ G R^^ \ {0}. (3.30) 

The next assertion is due to ^U\, [30] (see also [17] and, e.g., [M] Ch.6]). 

Proposition 8 The solution w G VJ/(S)^ of the problem p.29p admits the asymptotic form 

w{0 - (d(-V«)T$(e)^)Ta + [Di-V^miVVh + w{0. (3.31) 
|V|w(C)| <Cfcp-3-^ fcGNo = {0,l,2,...}, eGR3\V, (3.32) 

where a G R"^ and G while \a\ + \b\ < c\\g] L'^{dLu)\\. 

Remark 9 Formulae (|3.31|) - p.32[) can be derived from the integral representation of the solution w 
through the fundamental matrix (j3.30p . In this way decomposition p.32p is obtained from the Taylor 
formula in inverted variables ^|^|~^. Observe that the columns of matrices D{x)^ and d{x) in (|2.5[) and 
p.l3p form a basis in the linear subspacc of dimension 16 of columns linearly dependent on variables x = 
{xi,X2,X3). We emphasize that the matrix notation of elasticity relations combined with the polynomial 
property allow us to write the complete decomposition ofw{^) for |^| — ^ oo in a condensed and convenient 
form for further applications. 

Remark 10 Formula p.3ip contains the matrices d and D in p.lSp and ()2.5p . Let d^ {(^) , . . . , d"^ 
be columns of d{^) and let Di(^), . . . , Dq{^) be strings of D{^). Then we rewrite p.31|) in the form of 
strings 

j=i k=i 

Therefore, the asymptotic terms detached in p.3ip are but a linear combination of columns of the fun- 
damental matrix $(^) (with the coefficients ai, 02, 03 and 07; cf. (|3.13p ) and of the first- order derivatives 
of the columns (with the coefficients 04, 05, ag and &i, . . . , 6g ). 
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The columns d^,...,(f satisfy the homogeneous problem (|3.23p . However, the columns are not in 
the weighted space Vq'^(S)^ by the lack of their decay rate and, hence, {£,) are not solutions of the 
homogeneous {g = 0) problem (|3.29p in Proposition [51 According to the general method [29] such 
solutions are used to compute the coefficients in the asymptotic expansion (|3.31[) . We are going to 
use this method twice. First, we observe that the right-hand side g in (|3.23|) verifies the orthogonality 
conditions 

' d{O'^g{Ods^ = 0eR\ (3.33) 



Indeed, by (13.1511 . we get 

d{S.Vgms^ J = j D{n-md{Ods^ = -{A^e^)'^ J D{V^)d{OdC = 0. (3.34) 



Proposition 11 Under orthogonality condition (|3.33p . the column a G m p.3ip vanishes. 
The proof is commented in Remark 1161 

Let £ ^0^(0"' be a solution to the problem (|3.29p with the specific right-hand side 

g^iO = ^D{n-{OVA%,; (3.35) 

here j = 1, ... ,9, = ■ ■ ■ I'^j.g)^ is the unit column in R^, and 6j^k stands for the Kronecker 

symbol. Recalling the problem (j3.23p for the boundary layer term w^, we see that 

w\^)=W{Oe° (3.36) 

with the (4 x 9)-matrix function W composed from the columns W-^, . . . , of height 4, 

W = (W\...,W^). (3.37) 

By Proposition [TT] and the relation (|3.34p . we conclude the expansions 

W'iO'^ = ^ M,MW^)m^ + W^iO'^ (3.38) 
p=i 

where the remainders W^{£,) obey the estimates p.32p . The coefficients Mjp in p.38p form the matrix 
of size 9x9 

M^M{A^,uj) (3.39) 

which, in the analogy with |531 1381 1^ and others, is called the polarization matrix of the opening lo in 
the homogeneous piezoelectric space. 

As in Section 12.41 our study of general properties of (|3.39p relies on both formulations (|2.12p and 
(|2.22p of the piezoelectricity problem. Hence, we have to perform the same sign changes as in (|2.23p . 

M = ( ^^EM ^.^EE - = ) = ^^EM „^^EE ) • (3-40) 



Theorem 12 Entries of the modified polarization matrix M(^) satisfy the relation 

(M(=)),p = -Q'l_^iW^,WP; E) - (A°_)),pmes3t^, = 1, . . . , 9, (3.41) 
where Q^_'j is the quadratic form in (j2.22p with the matrix A^f^_^ = A(_)(0) (see (j2.23p and p.l6p ). 
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Proof. By ()3.35|) and (|3.15p . the sum yV^(C) = Dj{^)^ + verifies the homogeneous problem 

p.23p . In the method [55] these solutions play the same role as it was registered for the columns d^, . . . ,(f 
above Proposition [TT] We underline that the vector function 

>V('_) = {Wi'',Wt,Wi'',~W'^)'^ (3.42) 

verifies a homogeneous boundary value problem which is formally adjoint for (|3.23p and involves the 
differential operators L^-y^ and N^-y) constructed from and N'^ in (|3.17p and (|3.22p . respectively, by 
replacing with the transposed matrix (A")^. Clearly, L'^-j^{We) ~ L°(V^)* is the formally adjoint 
for the differential operator i"(V^). 

We insert and -j into the Green formula written for the truncated domain S/j = 5 n B/j and 
choose the radius of the ball B/j = {■f : |^| < R} such that the sphere Sr = 98^ envelopes the set uJ. We 
have 

{L'W\W^_^)s^ + (iV"W^^Wf_))a.us« = L^T)Wf_))s« + (W^^ iV(°T)Wf-))a-uSH. (3.43) 

Since L^YV^ = provides -^^-p-j W^_-) = 0, the integrals over in (|3.43p vanish. Furthermore, -^(V) {)'^\- 
0, ^Gduj. Thus, converts into 

{N°W^,Wl^)o^ = (W^^A^(°T)Wf_))s„ - (iV^W^-'', Wf_))s« (3.44) 

where N°{^,W^) = D{\^\^^^y A°D{0 on the sphere_§ij. 

Taking into account the estimates (|3.32p for and the concomitant estimates |V|VFP(^)| < 
Cpp^^"'', we obtain that the right-hand side I^fg/^^ of p.43p satisfies 

where means the asymptotic term detached in p.38p and Z)p(_)(^)^ is a column of the matrix D{(^)^ 
transformed according to (|3.35p . Understanding integrals over the ball B/j in the framework of the theory 
of distributions and using the Green formula, we obtain 

i^-f,Ht = (i"s^i?,V))»H-(^^^(T)^pV)W+o(^"') 

9 

= Y.^^^'iJ^p(-)(0D,i^iVm)dC + O{R'') (3.45) 

9 

= A/,,^9(-Vc)^p(-)(O^I?=o + 0(i?-i) 

9=1 



—Mjp for p = 1, . . . , 6, 

Mjp for p==7,8,9 

Here we have used that, first, £'p(-)(0 is linear in ^ and, therefore, LP^j-^D^^_^ = and, second, 

9 9 

L"(V,)I]-'(C) ■.= Y,M,,L\V,){D,{-V^)m''y ^Y.^'I,,D,{-V,y5{0 

q=l 9=1 

caused by the formula L"(V^)$(^) = 5(^)l4x4, i.e., by the definition of the fundamental matrix $. 

Let us process the left-hand side If^ft of (|3.43[) . Again integrating by parts, this time in the domains 
S and it follows that 

n!ft - (iv°w^^w^f_))oH-(ivX'^;(-))s- 

= Q\W\Wl_y^) + Q\D],Dl^_y,uj) (3.46) 
= Ql_^{W\W^- S) + (AO_)),p7nes3C., 
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where mes^uj is the volume of lu. Note that, first, the equality N^W^ = —N^Dj on duj is inherited 
from (|3.35p and (|3.15p . second, rf and —rf' imply the outward normals with respect to the sets S and 
w, respectively, and, third, 

Q\u,v^_)-^) = iA"D{V^)u,DiV^)v^^~,)s 

= iA°^^D{V^)u,DiV^)v)s^Q°^„^iu,v;E), (3.47) 
Q°^_^{Dj,D^;u;) = (A{'_)e„ e^)^ = (4„) ),pTOes3a;. 

Comparing (|3.45p and (|3.46p . we send R to +00 and obtain the desired relation (|3.41[) .B 
Theorem [T2] ensures the matrix in (|3.40p to be symmetric, in particular, M^^ = — (M^^)^. 

However, in contrast to the polarization matrix in elasticity (cf. [53l [38l |44] ) neither -/Vf(=), nor M enjoy 
the positivity/negativity property. In the case A^^ = Oexs the piezoelectricity problem decouples into 
the elasticity and electricity problems so that, 

< 0, M^^ > 0, = -(A/E^)^ = ©6x3, (3.48) 

provided, e.g., messU! > 0. We emphasize that in ()3.48|) M^^ is but the virtual mass tensor (see [48]). By 
the perturbation argument, the matrix Al has six negative and three positive eigenvalues, if the matrix 
A^^ is sufficiently small (cf. Section However, for arbitrary A^^, this property is still an open 

question. 

We have examined the first asymptotic term (j3.36p of the boundary layer type in the asymptotic 
ansatz p.l9p . By the representation (|3.38p (see Remark [T0|). we write the expansion of w^{^) for 
^ —>■ +00 in the matrix form as follows 

w\0 = iDiV^MCfVM^s^ + w\0. (3.49) 

The remainder obeys the estimates p.32p . 

Remark 13 Formula p.49p can be derived in the following way: 

p=l q=l 

\p=i 




3.4 The second term in the boundary layer 

The system of differential equations in S in the exterior problem p.23p for the boundary layer 
is homogeneous. This leads to relatively simple formulae (|3.36p and (|3.49p for w^. However, w'^ is 
determined from problem p.24p which enjoys the inhomegeneities both, in the boundary conditions and 
in the differential equations. Hence, the immediate objective becomes an inspection of the right-hand 
side —L'w^ for possible compensation and furthermore, an application of the same procedure as it is 
described in Section 13.31 However, the resulting asymptotic form (|3.59p of w'^ looks quite different 
compared to (|3.49p . 

By virtue of p.l6p and (|3.17p . the operator 

L'(C,V«)=i?(-V,)^ \J2^,AWiV^)] (3.50) 
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gets the following homogeneity property: 

L'it Vc)pV(^) = P^-'m, ? e \ {0}. (3.51) 

Here A S M, {p,0) are the spherical coordinates in M.^, p = |^| and 9 = p^^^ G Si, and f,ipG C^(Si)*. 
Thus, by means of ((3:36)) and ((338)) . ((3:32)) . ((330)) . we obtain that 

F'iO ^ ^L'i^,V^)w\0 ^ DiV^y ip-^^iO) + 0{p-'^), p^+oo, (3.52) 

while the formula can be differentiated under the standard convention V xO{p^^) = 0(p^'^^^). Due to 
the definition p.28p of the Kondratiev norm the right-hand side of (|3.52p gives rise to the continuous 
functional 



1/2 



\{F',v)^\ < c J p-MO\d^ <c(^J p-^dej 



\\p-\L\E)\\<C\\vXm. 



Thus, similarly to Proposition)?! we obtain the existence of a unique solution G Vq'^(E!)''' to the problem 
p.24() . Now, we need to examine the behavior of w-^{^) as p ^ +oo. According to [20] (see also (39l 
§3.5]), first of all, we have to determine the power-law solution 

ZiO=p-'Zi9) (3.53) 

to the system of differential equations 

i"(V^)Z(C) = p-^T{9) := D{V,yip~Hi9)), ^ e {0}, (3.54) 

with the right-hand side taken from ()3.52p . Note that, in general, the multiplier Z in ()3.53p may be 
linear in Inp but, owing to a special form of T, the next lemma proves the absence of the logarithm. 

Lemma 14 The system p.54p admits the power-law solution of form p.53p . whose angular part Z{9) 
is defined up to the linear combination ci^^ (9) + ■ ■ ■ + C4$^(0), where Cj G K are arbitrary and $^(0) is 
the trace on the unit sphere §i of the column $-'(^) in the fundamental matrix $. 

Proof. After separation of variables and rewriting the operator i^CVj) = p^^£,{9,We, pdp) in the 
spherical coordinates (p, 9) , the system p.54p takes the form 

£,{9,Vg,~l)Z{9) = T{9), 9 eSi. (3.55) 

By the Fredholm alternative, this system on the unit sphere has a solution if and only if the right-hand 
side J- is orthogonal to all solutions of the formally adjoint homogeneous system. Owing to (29) (see also 
[39) Lemma 3.5.9]), the formally adjoint operator for £(0, Ve,— 1) is nothing but £(t) (6*, Ve, 0), where 

p-''£^T){e,^e,pdp) - i{'T)(V?) = ^"(V^)*. (3.56) 

By virtue of the polynomial property ()3.27p . any power-law solution X(^) = p'^X{(^) of L^-j^{W^)X = 
in R'' \ {0} is a constant column in M'*. Thus, it suffices to verify the orthogonality condition 

J T{9)dse = e R-*. (3.57) 

Si 

Let i? > 7' > and let Q be the annulus : r < p < i?}. We have 

R 

^^{j^ j :F{0)dse= J p-^dp J T{9)dse = j p-^T{e)d£,= J D{W^f ip~^^{9))d^ 
= I D{p~^0^{p-'^{9))ds^~ I D{p-'O^{p-H{9))ds^^0. 
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We have used here the Gauss formula and the fact that the integrands aX p — R and p = r are equal to 
R^^D{9)^'${6) and r^^D{9)^"^{9), respectively, so that the integrals cancel each other. 

Thus, the compatibility condition (|3.57p holds true and the system (|3.55p admits a solution. It 
remains to recall that any power-law solution (|3.53p of the homogeneous system p.54p becomes a linear 
combination of the fundamental matrix columns. ■ 

To assure the uniqueness of the solution ()3.53p . we impose the condition 

J D{eyA°^{e, Ve, -l)Z{e)ds8 = O e M^ (3.58) 

Si 

where p~^'D{9,W0, pdp) is the matrix operator D{Wx) written, similarly to (|3.56p . in the spherical coor- 
dinates (p, 9). 

Now, we are in position to write an expansion at infinity for the second boundary layer term in ()3.19p . 



Proposition 15 The solution G Vq'^(^)'* of the problem (j3.24p admits the asymptotic form 

w\O^Z{i) + <^{£,)C + w\0, (3.59) 
\^'lw^{0\ < Cfe,^p-2-fc+/J^ fc e No, e e \ V, (3.60) 

where /J > is arbitrary, Z is a power-law solution of form (|3.53p and C £ US'* is determined as follows: 

C = -f{0)mes3Uj + J eR^, (3.61) 

J = [ i?(0)^^O^'^(Vc)(i?(Ve)$(O^)^ds^A/^e°. (3.62) 

Proof. The asymptotic expansion p.59p with a certain column C and the estimates p.60p result 
from [20] and [3^, respectively (see also [321 Ch.3]). We again employ the method proposed in f^Hj to 
evaluate the constant column C. Now, we use the Green formula in En for if^ and Cp = {Sp,i, ■ ■ ■ , ^p,4)^- 
Recalling (|3.24p . we have 



Iieft ■■=- 1 ejL' 


'w^d£_ - J e 






= J ejL°w^d( + 









. (3.63) 

= / ejN°w^ds^ =: I„ght- 

Here iV°(^, Vj) = D{9y A°D{V^) on the sphere with the imit normal vector 9 = p-'^i (cf. (|3?22)) 
and (|3.58p ). Similarly to the calculation (|3.45p . using (|3.59p and (|3.58p . we get 

Ir^ght = - / e], N^" Zds^ - J B^N^^ds^C + 0{R-^) = 

= j e^LHd^C + 0{R~^) ^Cp + 0{R-^). ^^'^^^ 

Br 

By integrating by parts, the last couple of integrals in lie ft turns into 

3 

ejL'w^d£.-J e^N'w^ds^ = J ej D{9)'^ J2^jAW{S/^)iD{S/^)^{0'^)'^dsiMe'^ + O{R-^). (3.65) 



Here we have applied the decomposition (|3.49p of together with the estimate p.32p for the remainder. 
Since its integrand is a positive homogeneous function in ^ of degree —2 (cf. p.30p ) the integral Jp over 
§i? in (|3.65p is independent of the radius R and becomes an entry of column p.62p . 
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The first couple of integrals in (|3.63|1 is equal to 

Here, the elementary formula (j3.18p has been taken into account. 

Now the limit passage R +00 in (|3.63p - (|3.65p furnishes (|3.6ip and (j3.62p .B 



Remark 16 Proposition can be proved by an application of the method ]29l in the same way as it 
is made in Proposition I J5I and Theorem We only mention that the columns , . . . , c?^ of the matrix 
d(^) in p.l3|) satisfy simultaneously the homogeneous problem p.23p and the formally adjoint boundary 
value problem in S with the operators L°-p^(Vj) and iVp-^ V^), respectively. 

3.5 The regular correction term 

Let us consider now the subsequent term in the asymptotic ansatz p.l9|) , namely the regular correction 
term m[x). 

By means of p.49p and (|3.59|1 . we have 



hw^{h-^x) + h^w^{h~^x) = h{S^{h~^x) + w^{h-~^x))+ 

+H\S^{h-^x) + w^{h-^x)) = h^S'^ix) + S\x)) + 0(/i4(|a-|-3 + |^|-2)) 



(3.66) 



where, according to p.30p and (|3.53p . we have set 

S'iO = {Di-S/^MO''yM^e°, = Z{0 + <i>(OC, (3.67) 

S'P(t^) =i-PS'P(e). 

Therefore, this is h^u{x) in the asymptotic ansatz (j3.19p that compensates the main part of a discrepancy 
produced by the boundary layer terms and w'^. 

Taking into account the equalities L'^S^ = and L'^S^ — —L'S^ designated in two last sections, we 
arrive at the following representation of the discrepancy in the system p.2p : 



t{x) = -L(x, V,)(x(x)(52(x) + S\x))) = 

= -[L,x]iS^ix) + S^ix)) ~ xix)iLix,V.,) - iO(V,) - L' {x,V ,))S^x)- (3.68) 
-x{xmx,V,)~L°iV,))SHx). 

Here [L, x] stands for the commutator of the differential operator L and the cut-of function x, i.e., 

[i, x\ = D{-V,)'^A{x)D{\7,x{x.)) - D{V,x{^))'^A{x)D{V,). (3.69) 

Recalling p.l6p and (|3.17p . in view of (|3.67p . we obtain that 

|f(a;)| < c|xr^ (3.70) 

We see that the regular correction term u must satisfy the piezoelectricity problem 

D{-V ^)'^ A{x)D{V ^)u{x) = f(.T), X e O, (3.71) 
D{n{x)yA{x)D{\/^)u{x) = 0, x G T^, u{x) = 0, a; G r„. (3.72) 

We emphasize that the sum hw^ {h~^ x) + h^w^{h~^x) in p.66p becomes of order h'^ only at a distance 
from the coordinate origin a; = 0. However, we have extended equations (j3.7ip over the whole domain Q 
because the singularity 0(|a;|~^) of the right-hand side f(x) is not too strong. In particular, by (|3.70p . 
the functional on the right-hand side in the integral identity 

Q(u,v;0) = (f,v)o, V e i?i(0;r„)4, (3.73) 
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serving for the problem (|3.7ip . p.72p (cf. (|2.12p ). is continuous due to the estimate 

1/2 / ^ 1/2 

\x\^\-i{x)\^dx\ I / N-|v(x)rdx| < 



i(f,v)oi<ciy' 



and the one-dimensional Hardy inequality mentioned above (|3.28p . Henee, in the analogy with Propo- 
sition [1] the Lax-Milgram lemma ensures the existence and uniqueness of the solution u e iJ^(ri; Fu)^. 
These observations complete the evaluation of all asymptotic terms detached in p.l9p . 

Remark 17 The singularity oft can lead to a logarithmical singularity of the solution u. However, we 
shall need only the following inequalities with arbitrary (3 > : 

|u(x)| < cp\x\-^, |V,u(x)| < cp\x\~^-^ (3.74) 

delivered by a result in \30\/ (see also \39\ %3.6]). 

For the further usage, it is convenient to rewrite the ansatz p.l9|l in a different form, namely 

u^{x) ^ u{x) + h^\5{x) + x{x){hw^{h-^x) + h^w^{h-^x)) + u''{x), (3.75) 

where, in accordance with (|3.66p and (|3.67p . 

V{x) = u(a;) + xix)iS^ix) + S^ix)). (3.76) 

In other words, we detach hS'^{h~^x) and S^{h^^x) from the boundary layer terms and attach them 
to the regular term u. Therefore, the remainder vP' in (j3.75[) stays the same as in the original ansatz 

Let us derive an almost explicit formula for (|3.76|) . To this end, let G{x, y) be the Green matrix for 
the piezoelectricity problem (|2.9p - (|2.1ip . i.e., 

D{-W^yA{x)D{W,)G{x,y) = 5{x - j/)l4x4, a; e 1), , . 

D{n{x)y A{x)D{V :,)G{x,y) = 0, x e F^, u{x) = 0, xe F„ ^ ' 

Of course, the relations p.77p are understood in the sense of distributions, so that, G S L^(r2)'*^^, 
G e L\dn)'^'"^ and 

{G,L^T)vh + iG,N^T)v)r^=viy), veC^{n;Tu)\ 

where the linear space C^{Q;Tu) consists of infinitely differentiable functions in Q which vanish on 
F„. Since A is a smooth matrix function inside of the ball Bjj (see p.7p ). the Green matrix is properly 
defined for y GBr (see [TTl [13]) and 

(x ^^ G(x, y) - $(a;, y)) e H\n)^''\ 

Moreover, G can be differentiated in the second argument and we set 

G"{x) = G{x, 0), G"(x-) ^ D{~Vy)G{x, y)\y=o- (3.78) 

By repeating the considerations in and around of Lemma 1141 we can detect that 

G° - $ e i^l(^7)4^^ G° - D{v.^)^ - z - K$ e H^in)^""^, (3.79) 

where K is a certain matrix of the size 9x4 with real entries and Z is such that Z{x) = Z{x)M^ e'^ (cf. 
p.50p and (|3.52p - (|3.54p ). Since, by definition of u and S'^, the vector function U verifies the boundary 
conditions (|3.72p and the homogeneous system (|3.7ip everywhere in fi, except at the point O. Let us 
now compare singularities in p.79p and p.76p to conclude that 

U(.t) = G°(a;)Af^£° - G^{x)f{0)mes3UJ. (3.80) 
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Remark 18 1. We emphasize that the differential operator Z?(— Vy) in p.78p is replaced by D(Vx) in 
p.79p . This is due to the evident relationship D(— Vj,)<i>(a; — y) ~ L)(Va;)$(a; — y). 

2. If A is a constant matrix, then the terms Z and K$ are absent in (j3.79p . in other words, their 
presence results from the variable coefficients of differential operator p.l7|) . Therefore, the column 
—f(0)mes3LU occurs on the right-hand side of p.80|) instead of the column p.61|) . To ensure that the 
additional column (|3.62p does not effect the form of the last term in (j3.80p . one may put = to see that 
then J = 0. A direct calculation leading to formula (j3.80|) can be found in 14 for the three-dimensional 
elasticity problem. 

Since the coordinate origin O is situated inside ujh, i.e., outside Qh (cf. Section [3.ip . the second term 
p.80[) in the new ansatz (|3.75p is smooth in the domain ^l{h), although the Green matrices (|3.78[) have 
singularities at O. 

4 Justification of asymptotics and analysis of shape functionals 

4.1 The justification of asymptotics 

The difference 

u'' ^u'' -u-xihw^ +h^w^)-h^u (4.1) 
(see p.l9p and p.75p ) satisfies the integral identity 

Qiu'',v;n{h))=T'\v), V e H\n{hy,Tu)\ (4.2) 

where is a certain functional (see, e.g., (|4.5p ). If the estimate 

\T'\>\ <ch''+''^^\\v;H\n{h))\\ (4.3) 

is proved, we could take v — w'' in order to conclude by using (j2.14|) that 

\\u'';H^{n{h))\\ Kch'^+^Z^ (4.4) 

In the sequel, it is shown, that the constants in (|4.3p and (|4.4p are independent of the small parameter 
h. 

To verify (|4.3p . first, we assume that v vanishes in the ball 82^/3, therefore, x^' ~ 0. Then, we have 
T''{v) = Q{u'' ~u- h^vL, v; n{h)) = Q{u'\ v; n{h)) - Q{u, v; n) - h^Q{u, v; f7). (4.5) 

Recalling (|3.6p . (|2.12p and (|3.73p . we observe that the support of the vector function (|3.68p satisfies 
suppi C 82^/3 (each term in (|3.68p contains either a cut-off function x supported in the ball, or its 
derivatives) and, hence, (|4.5p is null. 

Second, let suppv C \ ujh- We write 

F'^iv) = if,v)nih) - {ADiV^)u,DiV^)v)n,^,,-,^h^ADiV,)u,DiVMm 
-hiADiV,)ixw'), DiV,)v)n(H) ~ h\AD{V ^){xw^),D{^ M^(^h) (4-6) 
= : (/, v)a[h) - - h^l^ - hir ~ h^I^- 
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Since the vector functions u and u are smooth in B/f \ LOh, we integrate by parts and obtain 

3 

+ Y,{D{n^yx,A^s\ v)o^, + J\ (4.7) 
r = {D{n^)^AD{W,){u - D{xfe° ~ U),v)s^,+ 

3 

{D{n^)^ {A - A'' -Y,x,A^)e\v)e^, + {D{n^^y {A - A'')D{V ^)U,v)au.n. 

/" = (f,^^W)+?"> J'' = {D{n'yAD{V,)u,v)9^,. (4.8) 
To process the terms and we recah the inequality 

\x\~'^\v{x)\^dx < c\\v; H\n{h))f , (4.9) 



n{h) 



which is a consequence of the one-dimensional Hardy inequality (cf. [321 §4-5]) and the trace inequality 
(see [23]) 

\v{x)\'^ds^ < ch\\v; H\n{h))f , (4.10) 



where the constants c are independent of /i G (0, ho] and v. 
Now by (|4.10[) and (|3.74p . we readily derive that 

/i3|7"| < c/i^/i-i-'S / \v{x)\ds, < ch^-'^{mes2duJhy/^h^/^\\v;H\nih))\\ 

du^n (4.11) 
= Ch-P+y^\\v-H\VL{h))\\. 

Analogously, by means of (|3.16p . p. lip and (|4.10p . we have 

|7"| <c(/ii+" + /i2 + /i2) j \v{x)\ds, <ch''+^'^\\v;H^{n{h))\\. (4.12) 

We may choose /3 = 1 — a>0in order to equalize the final exponents of h in (|4.1ip and (|4.12p . 
Dealing with J™, we write 

- {AD{V^x)S\D{V^)v)nih)-{AD{V^)S\D{V^x)v)im 

+ {A^D{V,)w^D{V,){xv)hih)+I^, ^ ' 

= iADiV.,x)iw^~S^),DiV,)v)i^^^)-{ADiV^x)iw^-S^),DiV.,x)vhih) ^14^ 
+ {{A-A'^)DiV,)w^D{V,){xv))aw- ^ ' 

Here, we detach S^{h~^x) from w'^{h~^x) (cf. (|3.66p ) and commute twice the differential operator D{Wx) 
with the cut-off hmction x (see p.69p ). 

In view of (|3.16p and p.59p . the absolute value of the last expression in (|4.14p . multiplied by h'^ 
according to the definition of in (|4.6p . does not exceed the sum of the following two expressions: 

-3+13 



ch^ j \x\h-^(^^^ \D{V.){x{x)v{x))\dx< 
n\MRh 

/d^a,nn \l/2 (4.15) 

<ch^-P\ [ r^-^+^l^r^dr] \\v; H\n{h))\\ < ch^^Wv, H\n{h))\\ 



Rh 
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and 

ch^ J \x\\DiV,)w\h-'x)\\DiV,){xv)\dx 

1/2 

(4.16) 

<ch^Rh\ I h-^\D{\/^)wHO\''d^h'^ \ \\v;H\n{h))\\ 
< chy^\\v;H\n{h))\\. 

The radius R is chosen such that Bj? D uJ. Since the support of jDCVa^x)! belongs to the annulus 
^2R/3 \ B_R/3 where, according to (|3.59p . 

\w\h-^x) - S\h-^x)\ + \V,iw^{h-^x) ~ S^{h-^x))\ < ch^-f^, 

the remaining terms in (|4.14[) . again after multiphcation by h^, are bounded by ch'^^'^\\v; H^{^1)\\ while 
we may set /? = 1/2 to achieve the same exponent as in (|4.15p . In other words, for p ~ 2, we now have 

hp\i;;'\<chy^v;H\n{h))\\. (4.17) 

By formulae (|3.49p , (|3.32p and (I3.66P , the similar argument leads to the estimate (|4.17p for the remainder 
in the representation 

(A"i?(v,V\i?(v,)(x«))o(h) + E'=i(2^j-^''^(v.)w\i?(v,)(x«))o(h) + /r. ^ ■ ' 

Now, we arc in position to conclude the estimate (|4.3p for the functional F'^ in (|4.2p , (|4.5p and (|4.6p . 
To this end, wc list several facts. First, the inner product (/, on the right hand-side of (|4.6p 

cancels the same product in (|4.7p . Second, the equality 



is inherited from the definitions ()3.68p and ()3.69p . Third, we make the coordinate dilation .t i-^ ^ = hr^x 
in the first couples of terms on the right hand-side of (|4.13p and (|4.18p . simultaneously multiplying the 
terms by and h, respectively. Noting that S''^{h~^x) = hPS'''{x), p — 1,2, we sec that these couples 
and /i'^(f , v)j7(-;j) annihilate. Finally, we recall the integral identities (|3.29p . serving for the problems 
p.23p and (|3.24p . and after the substitutions .t i-^ ^ and i— *■ x(^C)''^(^C)i we detect all terms in the 
identities on the right hand-sides of ((i?7)l . (|iT^ and ((i7T5)l . Thus, 

Ff^iv) = r + h^J"^ + hi]" + h^i^ 

and the inequality ([43]) holds true by virtue of ([4T2l) . ([4TT|) and (|4T7l) with p = 1, 2. We notice that 
the lowest exponent a + 5/2 of h occurs in (|4.12p because a G (1/2, 1) and a + 5/2 S (3, 7/2). 
We now formulate the result. 

Theorem 19 Let all assumptions in Section \3.1\ be valid, in particular, the inclusion (|3.9p with a S 
(1/2,1). Then the solution of the piezoelectricity problem (|3.2p - (|3.5p and its approximation con- 
structed in Section\^ are in the relationship 

\\u^ -u-h^n^ X{hw^ + h^w^)-H^{n{h))\\ < ch°'+^/^N, (4.19) 

where the constant c is independent of the parameter h G (0, /iq] and the right-hand sides f , g while 

N=\\f- L\m\ + L\dm\ + II/; C2'"(Bfl.)ll- (4.20) 

In asymptotic estimate (|4.19p u stands for a solution of piezoelectricity mixed boundary value problem 
(|2.9p - (|2.1ip in the intact body and u for the main regular corrector which is a solution of problem 
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p.7ip . p.72p in 51 and admits the representation (|3.76p with the smooth addendum U and the singular 
components (j3.67p . The boundary layer terms and uP' , which are as well present in (j4.19p . are given 
by solutions to the exterior problems (j3.23p and (|3.24[) . decay at infinity and take asymptotic forms (j3.49p 
and p.59p . respectively. 

Remark 20 The obtained estimate ()4.19p is asymptotically sharp, in particular, it satisfies the "first 
omitted term" rule. Indeed, for the smooth data A and f , the subsequent asymptotic term in the ansatz 
p.l9p is h'^xi^)'^^if^^^^)! ihe H^{Q{h))-norm of the latter term is just 0{h''P). This bound appears 
in (|4.20p if a 1 — 0. Moreover, the estimate (|4.19p holds true when the last addendum in (|4.20p is 
changed for ||/; C^'"^(Bfl)|| with any ai G (0,1). If the right-hand side f € C^'"(Bfl)'* in the equations 
p.2p is not sufficiently smooth, e.g., 

fix) = /"(x) + \x\'+"f\0), f G C-(B^)^ /I € C-(§i)4, 

then the asymptotic ansatz p.l9p gains the boundary layer term /i^+"x(2;)u'^+"(/i^^a;) with the Sobolev 
norm in fl{h) of the some order /i"+5/2 qji ^/jg right hand-side of (j4.19p . 

A direct calculation show that 

h^\\xw^-H\^{h))\\ = 0{h^+'/% J = 1,2, (4.21) 

and, therefore, in view of the relation q + 5/2 > 3 (see p.9p ). the i7^(Jl(/i))-norm of each of the detached 
asymptotic terms in (|4.19p (cf. (j3.19p and ()3.75p ) is of order h" with s < 3. In other words, Theorem [T9l 
justifies the constructed asymptotics of solution it'', indeed. 



4.2 The energy and the electric enthalpy 

We proceed with energy functional (|2.17p . assuming for simplicity that the volume forces and the volume 
charges are absent, i.e., / = on the right hand-sides of (|2.9p and (|3.2p . Then, integrating by parts and 
taking into account formulae (|3.75p and (|4.4p . we have 

Uiuf^; n{h)) = \{D{nyAD{^^)u\ u'^)r^ - [g, u^)r^ = -i(.g, u^^)t. 

= -\{9.u)v. ~ \h^{g,v.\^^0(h'^^^l^) . ^^■''> 

Let e H'^{n]Tu)'^ and 0^ G ifi(fl;r„)'^ imply the solutions of the problem ([ ^ - (pUj) with the 
right-hand sides 

0^ = (.9r,52^,52'',O)T, 0E^(O,O,O,gf)^. (4.23) 
Using the representation (|3.80p with /(O) = and the modified column U(_) (see (|3.42l) ). we obtain 

(fl^,U)r„ = (0^,U(_))r„ = (i?(n)TAi?(V,)o^,U(_))r„ 
= (o^,i?(-V,)T^Ti?(V,)U(_))r^ 
= {s'yM{^^,D{-\/,yA-^D{\/,)Gl^)n 

= (eO)TA/(o^,(i?(V,)Tj(x))(_))o = -(eTA^e^_), ^ ' 

(0E,U)r„ = (fl^U(„))r„ =-(oE,Z?(-V,)TATi?(V,)U(_))r„ 
= (£0)^Mef_), 

where e'^ = D(V^)t)'^(0) and c^^ = D(y^)v^{0). Here, we apply formula ([XT?]) for the derivatives G° of 
the Green matrix G in (|3.77p . We emphasize that 

(/,GO_))o + (.g,GO_))r^ =-(Z?(V,)«)(_)(0) = -e5'_) (4.25) 

because entries of G*^ are given by the derivatives of columns of the fundamental matrix G{x,y) with 
respect to the second argument, and satisfies the problem 

Di-V,y Aixy D(y,)G^^){x, y) = 5{x ~ 2/)I(_), x e n, 

D{n{x))'^ A{xy D(y,)G^^){x,y) = 0, x € T,, G(_)(.t,2;) = 0, x G r„, 
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where — diag{l, 1, 1, —1} (cf. problem (j3.77p ). 

By (|4.22p and (|4.24p . the following representation is valid: 



U{u''; nih)) - U{u- n) = _(£0)TM(e^_) - cf_^) + 0{h''+^'^). (4.26) 
At the first sight, ()4.26p looks like (j2.ip . however this impression is wrong. 

Remark 21 The decomposition u = + is only a mathematical device in our analysis, since 
in a smart material it is impossible to distinguish between the strain columns = D{'S/x)v^{0) and 
= D{Wx)'O^{0) generated at the point O by the external mechanical loading and the electrical 
surface charge in (|4.23[) . Surely, one can measure only the sum = DCS/ x)u{0) resulting from 
complete external action and standing as the first term on the right hand-side of (j4.26p . This unususal 
property of the energy functional should always he taken into account and ignoring the above observation 
on the decomposition u = 0^ + may provoke for misleading physical conclusions. Similar situation 
occurs for example for a crack in a piezoelectric medium. Applying the Griffith energy fracture criterion, 
in 14-6}/ the energy release rate at the crack tip is expressed in terms of stress intensity factors, i.e., local 
characteristics of the elastic/ electric state at the tip. In \2S^ a mistake in a calculation (formulas (33.23) 
and (34.4s) in 146] pages 296 and 312]; cf. comments in \22^ } was discovered and a non-local formula for 
the energy release rate of type (j4.26p was derived rigourously and justified. The non-local character of the 
energy release rate means that the energy functional U {u] fl) cannot be employed for a fracture criterion 
and in the Griffith criterion must involve the electric enthalpy ( cf. U5l \46^ for an interpretation from 
the view point of solid fracture mechanics) . 

The difference 

e - e = (0 , Gr )r„ - (0 , )r„ 

ought to be regarded as a global characteristics of the mechanical electric state of the body fl and, 
therefore, formula (|4.26p has a different physical meaning compared to (|2.ip and (|4.33p below. 

Let us now compute the increment £{u^\ ^l{h)) — £{u; f2) of the mechanical enthalpy determined in 
(|2.20p and (|2.2ip . Returning back to the general case f 0, we obtain 

Eiu'^^nih)) = ^{AD{Vx)u^DiVx)u1_^h(h) - (/,«f_))n - (.9,<_))n 
= ^iDi-VxVAD{Vx)u^u'(_^)nih) + ^{Diny AD{Vx)u'\u>(_^)r^ (4.27) 
- i9,u'l_^)n = -i(/,wC_))o - i(5,<_))r„. 

As above, we have 

(5, u\_^)v. = (5, «(-))r„ + h^g, U(_))r„ + 0(/j"+^/2). (4.28) 
Furthermore, in view of representation p.75p we derive 

{f,u'(_^h - {f,u^-)h + h'{f,Vi-)h + 0(/i"+'/'). (4.29) 
according to inequality (|4.4p and the following relations 




LUh 



h\{f,xw')nih)\<ch j (l + ^)-Vdr<c/i4|ln/i|<c/i"+5/2, (4.30) 


h\f,xw^)niH)\<ch' j (l + ^)-2+Vdr<c/.4-^<c/."+^/2. (4.31) 




23 



In the estimation (j4.30p wc have used the formulae (|3.32p and (|3.59p for and w'^ together with the 
demanded inclusions a E (1/2, 1) and S E (0, 1/2). 

Now formulae ((4281), ^j^^ and ((3?80)) . (|4:^ convert ([4271 into the form 

-lh^{e'^yMiif,Gl^h + ig,Gl^)r^) + 0{h-+'/^) = (4.32) 

- W-(0)^/(0)mes3C.,) + i/i3(eO)T^^^.o_^ ^ 0(/j"+5/2) 
= i/l3(eO)TM(^)£0 + Ci(/i"+5/2). 

Here, we have taken into account that, first, Me^_-^ = M(^)e° according to the definition of in 

p.40p and, second, ii_(0)^/(0)mes3(jj/i = (/, U(_))c^^ +0{h^^°') due to the smoothness properties (|3.10p 
and (|3.1ip of / and u. 

Let us formulate the result obtained in (|4.32p . 

Theorem 22 The electrical enthalpy (j2.20p admits the asymptotic expansion 

f (w''; r!(;i)) = £(«; n) + i/i3(eO)TM(=)eO + 0(/i"+^/2), (4.33) 

where and u imply solutions of the piezoelectricity problems p.2p - (l3.5p and (|2.9p - (|2.1ip . respectively, 
= Z)(Va;)u(0) is t/ie strain column (j3.12p and M^^j = Af^^-) (A*^, w) is f/ie modified polarization matrix 
which is a symmetric matrix of size 9x9 /see formulae (j3.39p . p.40p and r/teorem lig)) . 

Note that in contrast to the energy functional (|2.16p the electrical enthalpy has the topological 
derivative 

i(i?(V,)w(0))^M(=)(A",u.,0^(V.)w(0) (4.34) 

expressed in terms of local characteristics of the elastic/electric state in the entire body and of the 
shape of the small void Uh- Owing to representation (|3.4ip . we emphasize that the polarization matrix 
p.39p enjoys the homogeneity property M{A^;ujh) — h^M{A^\uj) which has been used in the passage 
from Km to ((4:341) . 

Notice, that exponent 3 in the factor is conform with formula (j3.4ip for polarization matrix which 
contains the volume messw of the void. 

4.3 Shape functionals and the adjoint state 

Possible applications of asymptotic analysis performed in the paper include inverse problems, optimum 
design and shape optimization. We refer the reader to [HIHIISI for numerical results of shape and topology 
optimization by an application of the levelset method. In the levelset method the topological derivative 
of a specific shape functional is employed to detect the regions of the hold-all domain to include voids in 
order to improve the value of the functional to be optimized. The numerical method turn out to be very 
efficient in two spatial dimensions compared to the pure levelset strategy. Another application with, it 
seems, very high potential for numerical solution are all types of inverse problems to detect imperfections 
within a geometrical domain on the basis of boundary observations. However, in inverse problems it is 
required that the data imply the unique identification of the imperfection. This property is unknown, 
in general, for the strategy which is based on the asymptotic analysis in singularly perturbed domains. 
In particular, it is an open problem how to identify an imperfection from the observation of a finite 
number of cigenmodes (eigenvalues), which seems to be a natural and efficient way to solve the problem. 
The difficulty of such an approach is hidden in the fine properties of topological derivatives which are 
still to be investigated, for example that some sufficiently large set of observations in the mathematical 
model based on the asymptotic analysis leads to the uniqueness of the position of imperfection. It means 
that the derivation of topological derivatives is far from being sufficient for the practical applications of 



24 



the promising tool of shape and topology optimization and identification. We can consider below some 
specific examples of shape functionals, the other possibilities, including the spectral problem require 
some additional work to derive the asymptotic formulae. 

We return to the analysis. Recalling the Sobolev embedding theorem iJ^(ri) C L^{fl) in M.^, we 
assume that the density J in the shape functional 



J{u;n) = J J{u{x);x)dx (4.35) 
n 

satisfies the following restrictions: 

|J(a;x)| <c(l + |a|*), (4.36) 

\J{b; x) - J{a; x) - J'{a; x)^(6 -~a)\<c\a~ b\'^{l + \a\'~^ + |6|*~^), (4.37) 

\J{b;x) - J(6;0)| < c\x\''{l + |6|*) (4.38) 
where x G ft, a and b are arbitrary columns in W^, and the vector function J' is subject to the conditions 

|J'(a;x)| <c(l + |a|*-i), (4.39) 

I J'(a; x) - J'ib; y)\ < c{\a - &r(|a|*-^ + l^l*'^) + \x ~ ypi\a\' + \b\% (4.40) 

while 

i G [2,6), 7G (0,1). (4.41) 

In other words, along with the restrictions on the growth of J and J', the integrand J is differentiable 
with respect to the first variable and Holder continuous with respect to the second variable. Moreover, 
J' is Holder continuous in both arguments. Inequality (|4.36[) ensure that functional (|4.35p is defined for 

Remark 23 Simple examples 

R{x)\u'\x)-u{x)\^dx and J R{x)\u'' {x)\^ dx 
n{h) n{h) 

are related to the least square method and satisfy the above requirements with t = 2 and 7=1 for 
R € C^'"{fl). In addition, for g = m the boundary conditions (j2.10p and p.3p . the work of external 
forces 

' f^{x)'^u''^{x)dx+ f f^{xyu''^{x)dx 



o(/i) a(/t) 
and the electric enthalpy ( cf. (14.27^ ) 

£(w^^^(/l)) = -i J /^(x)^«''^(x)dx + i J fHxfu'^Hx)dx 

n{h) Q{h) 

readily display another examples. 
We consider the difference 

J{u'';n{h))-J{u;n)= J {J{u^{x);x) - J{u{x);x))dx + j J{u{x)-x)dx. (4.42) 

n{h) 

and, owing to (|4.37p and (|3.75[) . obtain the formula 

\J{u^^{x)-x) - J{u{x)-x) - J'{u{x)-xy{h^\J{x)+x{x)Y.]=ih'^^{-h)+u\x))\ < 

<c(/i6|U(x)|2+x(x)2V'^^/l2j|iCJ(i)|2 + |il'^(^)|2)(l+|y;»(.^)|t-2^|y(^)|*-2)_ ^ ' ) 
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Recalling the estimates ()4.4p , (|4.39|) and applying the Holder inequality with the index couples {p, q) 
(5/6,6) and {p,q) — (3,2/3), we obtain 



rH.);.r='-w<i.<c/(i + M.,ftis'.wi..< 



n(h) 



n{h} 



<c{l + \\u;L'^{n)\\''')\\u^;L%n{h))\\ < c\\u''; H\n)\\ < c/i"+5/2, 



\u''Wl + \u"r' + \ur')dx <c J \u''\'{l + \u''\^ + \u\'')dx< 
n{h) n{h) 

< c\\u''; L^{n{h))\Wl + L%n{h))\\^ + \\u; L%n)\\^) < ch'^°'+^ 



Similarly, 



/i^ J \\J{x)\^{l + \u''{x)\*-^ + \u{x)\'-^)dj 



< ch". 



However, because of the singularity |U(a::)| = 0{\x\ '^), we use in the ball Mr/ the Holder inequality with 
the couple 

6 6 



(4.44) 



to derive that 



|U|^(1 + + |it|*-2)d.T < c/iM / r-^r'dr 



x{l + \\u'';H\n{h))\\*-'^ + ||M;i7i(17)||*-2) < c/i^-t/s. 



Wc deal with the boundary layers in the same way as in (|4.15[) and (|4.16p . Outside the ball we 
apply the inequalities p.32p and (|3.59p even much rougher ones, to conclude by the Holder inequality 
with the index couple (|4.44p that 



xi^W (l) f (1 + \u''{x)\^ + \uix)\^)dx 



< 



< ch^ 



(4.45) 



Inside the ball lS>B.h the Holder inequality gives 



X 

w' I — 



{l + \u''{x)\' + \u{x)\')dx 



Note that =37 < 6 due to (|4.4ip and, therefore. 



\\w^;L^^{Mr\uj)\\ < c\\w^;H\Mr\uj)\\ < c\\w^ ; VoHE)\\ . 
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Although, the faster rates of decay of the remainders and w'^ (cf. (|3.66p ) are not used in the 
estimation (|4.45p , the rate of decay becomes an important ingredient of the inequahties 

h^ J J'{u{x)-xyx{^)w'({)dx\<ch''^ 3 = 1,2, 
n(h) 

its derivation is much simpler, though. A simplification originates from the relation \ J'(u{x); x)\ < const 
for X G suppx C Bfl" so that one may repeat the calculation (|4.30p . 
Finally, we write 

Rh 

h^\ [ J'{u{x);xyV{x)dx\ < ch^ [ r-^r^dr < ch'^ 



and, in view of ([XTTj) and (gHOl), 



J{u{x); x)dx — h J{u{0); 0)TOes3W 



< c/j3+™™{a,7} 



Everything is prepared to derive a formula of type (|2.ip for the shape functional (|4.35p . 
Theorem 24 Let the assumption formulated above hold true. Then the asymptotic formula 
J{u'';n{h)) = J(u;n) + h^{{J{u{0);Q)-P{0yf{0))mes3Uj 



(4.46) 



is valid where P G H^{^;T^i)'^ H C^'™*"'f"''''^(]B^/)^ is a solution of the formally adjoint piezoelectricity 
problem 

D{-v.,yA{xyD{v,)p{x) = j'{u{xy,x), xen, 

D{n{x)yA{xyD{\/^)P{x) = 0, x £ r„,P{x) = 0, x e T^. ^ '> 

Proof. The calculations performed above provide the relation 

h-\j{u'';n{h)) - J{u- n)) = J(u(0); 0)mes3C^ + ( J'(m), U)o + 0(;i™™{7."-i/2,3-t/2})^ 

We recall the representation (|3.80p where is the Green matrix, i.e., a solution to the problem (|3.77p . 
The Green matrix and its derivatives help to calculate the solution P of the formally adjoint problem 
(|4.47[) and the derivatives G° (see p.78p ) deliver the column DiV x)P{x) at the point x = 0. In other 
words, we write 

(J'(w),U)a = {D{-V,YA^D{V^)P,G^)nMe'' 

- mesMD{-^xYA^D(^,)P,G%,f{Q) 
= iP,DiV,ydMe°)n-mes3ioiP,6f{0))a 
= -{D{y^)P(Q)Y Me^ - P(0)^/(0)mes3W. 

We again used the Dirac mass 5 in the framework of the theory of distributions to compute the expression 
(|05D . 

Finally, in order to justify our calculations we make the following comments. By assumptions (|4.36p . 
(|i:iT|) and (11301), (EH), the functional 

H\n;Tuf3v^{J'{u),v)n 
is continuous and J'{u) e C°'™™'f"''''''(IBi?') with any R' < R. Thus, the same arguments as in Sections 



and 13.11 guarantee the existence of a solution P to the problem (|4.47p which is twice differentiablc in 
the vicinity of the point a; = 0. These observations make all calculations justified. ■ 
The topological derivative of the functional J', i.e., 

T{u, to) = ( J(m(0); 0) - P(0)T/(0))mes3t^ - (i?(V,)P(0))^Me°, 

is non-local since it involves the adjoint state P in (|4.47p which depends on the solution u of the 
piezoelectricity problem in the entire domain Q. 
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4.4 Example 



It turns out, that the so-called weak interaction is quite common feature of piezoelectric materials. 
Therefore, we are going to present an example for such materials. Assume that there is a weak interaction 
between the mechanical and electric fields. This means that in the decomposition 



A - A(o) + , (4.49) 



06X3 A . ( ©6x6 -^"^ 



"^("^ - \ ©3X6 ) ' ^(^) - V A^^ ^3X3 

the entries of matrix are much smaller compared to non trivial entries of the matrix A(o) • It implies 
that in the first order approximation the piezoelectricity problem is decoupled into two problems, the pure 
elasticity problem with the stiffness matrix A^^, and the pure electricity problem with the permeability 
matrix A^^. 

We are going to evaluate the main correction terms in the asymptotic expansions of characteristics 
for the piezoelectric bodies fl, S and (see Sections 12.31 13.31 and l3.ip . We point out that to evaluate 
the main asymptotic terms in all formulae given below, it is sufficient to solve only the pure elasticity 
and the pure electricity problems. 

Remark 25 Since we always deal with the first order asymptotic corrections, the introduction of any 
small amplitude parameter neither makes formulae more transparent, nor contribute to the exactness 
of presentation. We emphasize that, in contrast to the preceding sections, the perturbations here are of 
regular type, which means that the justification of obtained formulae relies upon the standard argument 
of convergence of Neumann series. In order to simplify the notation, in the sequel the second order terms 
are always denoted by dots, starting from (|4.50p . 

We proceed with the solution 

u(a;) U(o) + U(i)(a;) + . . . (4.50) 

of the problem (|2.9p - (|2.1ip . In view of (|4.49p . the displacement vector and the electric vector u^^^ 
verify the problems 

i?^(-V,)^A^^(x)i?^(V,)zi^o)(x-) = /^(x), xen, (4.51) 
D^(n(a;))TA^^(x)i?^(V,)7/^o)(a;) - .9^(x), x e T,, ^^^^(x) = 0, x e r„, 

-VjAE^V,Mfo)(a;) = f^{x), x e n, (4.52) 
A^^ufg-^{x) = g^{x), x e T„, uf^-^ix) =0, x S r„ , 

and can be determined separately. Inserting (|4.50p and (|4.49p into (|2.9p - (|2.1ip . we arrive at the problem 



D(-V^)^A(o)(a;)I?(V^)M(i)(a;) = D(-V,0^(i)(a;)i?(V:,)u(o)(a;), x € Q, 
I?(n(x-))^A(o)(a;)I?(V,)?/(i)(x-) = i?(n(.T))^ (.T)i?(V,)ii(o) (x), x & T,, (4.53) 

U(i)(a;) = 0, X eVu . 

This problem is decoupled as well, however, its solution manifests the interaction between electric and 
mechanical fields, since the displacement vector m^'j depends only on the main part w^qj of the electric 
potential and, in the same manner, m^j^^ depends on u^y 

In order to complete the asymptotic formulae, in the same way as in the previous sections, we also 
need the expansion for the polarization matrix 



28 



M = M(o) + M(i) + . . . , (4.54) 

Wc emphasize that the matrices M(o) and inherit the block diagonal structure of A(o) ^^nd 

the block-anti-diagonal of ^(i), respectively. The same structures are kept by all matrix objects, in 
particular, the fundamental matrix takes the form 

(4.55) 




Here, is the fundamental matrix for the elasticity matrix operator D'^ {-\/ A^'^ D'^ {\/ ^) and <i>^pj 
is the fundamental matrix for the scalar operator — VjA^^V^. Furthermore, Af(o) -^(o) ^^'^ 
virtual mass matrix and the elasticity polarization matrix for the cavity a; C M^, which are negative 
definite (see [48] and [H [44]). 

It is convenient to proceed with the matrix solution (|3.37p which, according to (|4.49p and (|3.38p . 
enjoys the expansion 



with 




(4.56) 



©3x3 \ 03X6 W^(7f 



W{0 = {MD{S/^)m'^V + O(l^r') = (4.57) 
(M(o)i^(Ve)<i>(o)(e)^)^ + (M(o)i^(V^)$(i)(e)T + A/(i)7^(Ve<i>(o)(e)^)^ + • • • + 0^-') . 

The correction term in (|4.55[) is a power-law solution of form (|3.53p for the system of differential 
equations 

i?(-V^)^A|'o)i?(V4)$(i)(C) = Z?(Ve)^4)i?(V^)$(o)(0 , e e M3 \ {0} , (4.58) 

(cf. (|3.54p ). By a general result in [20] (see also [39l Lemmas 3.3.1 and 3.5.11]), the solution can 
depend linearly on ln|^|, however, the same argument as in the proof of Lemma 1141 ensures that <i>(]^) 
is positive homogeneous of degree -1 according to p.30p . The solution <&(i), which is defined up to the 

'(o)( 

fixed such that 



linear combination <&(o)C' of the fundamental matrix columns with the constant column C G R^, can be 



J i?(V^)^A°i:i(V5)$(i)(Ods5 = e . (4.59) 

The exterior problem for the correction term in (j4.56p takes the form 

i?(-Vj)TA^„)i^(Ve)W^(i)(C) = i^(Ve)TA|'i)i?(V^)W^(o)(e) ,e e 5 , (4.60) 

i^(n-(0)TAOo)Z?(V4)M^(i)(0 = -Din'^iOVAl^DiV^)W^o)iO ,(&du;. (4.61) 
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Since, owing to (|4.57p . we have 

W^o)iO = (M(o)^(Vc)$(o)(0^)^ + 0(|^|-3) , 
the right-hand side i^(i)(C) in (j4.60|l admits the decomposition 

i^(i)(0 = i?(V,.)^4)i?(V,.)(i?(Ve)$(o)(C)^)^A/({,) = 

3 ^ 

= E^^(^«)^4)^(^?)*(0)(^)^^K)^A^J) = (4.62) 



9 

3 



with the remainder = In (|4.62p . the equation (|4.58p has been apphed. Comparing 

with (|i37)) . we set 

WiDiO = + (M(o)i?(V^)<i>(o)(C)^)^ • (4.63) 

Recall that lo contains the origin ^ = 0, therefore, the last term in (j4.63p is smooth in S. As a result, a 
new exterior problem is obtained, with the right-hand side which decays sufficiently fast at infinity, 

Ve)^ AOo)i?(V5 )W?(i) (0 = F(i) (0 , ? e 5 , (4.64) 
Z?(n-(C))^AOo)i^(V4)t?(i)(0 = G(i)(0 ,C e au; , 

where 

G(i)(e) = D{n'^iOVAl^D{V^)W^om- (4.65) 
-I?K(e))^4o)i?(V«)(A/(o)I?(V^)$(i)(e)^)T . 

Now, the decay of G(i)(^) can be used, indeed, by Proposition [5] (see [50] and [321 Theorem 3.5.6]) and 
the calculations (|3.34p . p.33p . the solution e Vq'(S)"' admits the asymptotic form 

- (M(i)i?(Ve)<f>(o)(e)^)^ + iy(i)(0 , (4.66) 

where the remainder is subject to the estimates (|3.32p with the majorants CkP^^^'^^^ ((5 > is 

arbitrary) and the notation used for the derivatives of the fundamental matrix <I>(o) is matched with 
formulae fTFT)) and (g^. 

In order to evaluate the correction term in the expansion of the polarization matrix the method 
|29j is employed, here we recall that the columns of the matrix 

>V(o)(-)(0 = Di-m'^ + W(o)i-){0 (4.67) 

(cf. (|3.4ip ) are formal solutions to the homogeneous problem (|4.64p . By the Green formula in S n B^, 
we obtain 



W(o)(_)(C)^F(i)(e)de + J W(o)(-)(C)^G(i)(C)ds4 = (4.68) 

(i?(i)(o^i^(ier'o^4o)(-)^(e)W(o)(-)-wj)(_)i?(ier'o^^(o)(-)^(o%)(0)ds^ 

+0(i?-i)--M(i)(=)+0(i?-i) . 
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Wc have here repeated the computation (j3.45p based on the representations (|4.66p and (|4.67|) . The 
integrand on the left-hand side of (|4.68p is of order |^|^^ and, hence, the integral over S converges and 
the formula 

CO Ad^^ \ P ~ r ~ 

Jem ^^^^^ ]=-J W(o)(-)(e)^^^(i)(OrfC + J W(o)(-)(e)^G(i)(Od5« (4.69) 

together with (|4.62|l - ()4.64p expresses the matrix M(i) (cf. the definition (j3.40p ) in terms of the matrix 
^ME ^ (^EM)-i g^j^j ^Yic special solutions W^, . . . and W'' , W^, of the pure elasticity and the 
pure electricity exterior problems in S. Theorem [T^ shows that (M^^)'^ = —71/™. 

The formulae derived above can be used, e.g., to obtain the topological derivative of the electric 
enthalpy 

T£{u;uj)= (4.70) 
= i/i3((i?^(V.)^^o)(0))^A/(^)i?^(V.)z.^o)(0) - V.4)(0)^A/fo)V.«fo)(0)) + 

+/i3((i5^(V.)n^o)(0))^A/(o)^''(V.O<i)(0) - V,,4)(0)^Affo)V.4)(0))+ 
+/i3v,^.fo)(0))^M™i?^(V,)^.^i)(0) + . . . , 
where and A^^g) elasticity polarization matrix and the virtual mass matrix for the cavity 

while Mf^ = - (^MlffY is expressed in glHl)- 

Even the main term (with the factor ^h^) of the topological derivative (|4.70p has no sign, that is, 
in contrast to the forms of topological derivatives of the energy functionals for the pure elasticity and 
the pure electricity problems. The correction term (with factor h^) in (j4.70|) depends on two specific 
ingredients, namely, the correction term M^^ in polarization matrix (see (|4.54p and (|4.69p ). and the 
correction terms u^j, u^-^^ for the combined mechanical and electric fields. 

Remark 26 All the attributes in the above formulae can be given explicitly for some canonical shapes, 
including balls, ellipsoids and elliptic cracks in three spatial dimensions, and some other shapes in two 
spatial dimensions (see 148}/ and \5S\. \3'A \25[ and others). 

Remark 27 The case of g^ = QJ^ = has a very clear physical meaning (i.e. one gets an electric 
sparkle when pressing the lighter button). Then, in notation of Section \4-.2[ 

u^=u, u^ = 0, e^=e°, e^ = 0, 

thus, by relation (j3.40p . we can conclude that the topological derivatives in (j3.40p and p.46p of the energy 
and electric enthalpy functionals coincides one with another. In general, this identity is false, and can 
be misleading for the choice of governing Gibbs' functional for piezoelectric body (cf. Remark\21\). The 
relations between the topological derivatives for elasticity and piezoelectricity are easy to established, since 
the topological derivative for piezoelectricity can be viewed as the difference of that for elasticity and of 
the other for electricity. 
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